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In order to assess future damage caused by natural disasters, it is desirable to estimate the damage 
caused by single events. So called damage functions provide - for a natural disaster of certain 
magnitude - a specific damage value. However, in general, the functional form of such damage 
functions is unknown. We study the distributions of recorded flood damages on extended scales 
and deduce which damage functions lead to such distributions when the floods obey Generalized 
Extreme Value statistics and follow Generalized Pareto distributions. Based on the finding of broad 
damage distributions we investigate two possible functional forms to characterize the data. In the 
case of Gumbel distributed extreme events, (i) a power-law distribution density with an exponent 
close to 2 (Zipf’s law) implies an exponential damage function; (ii) stretched exponential distribution 
densities imply power-law damage functions. In the case of Weibull (Frechet) distributed extreme 
events we find correspondingly steeper (less steep) damage functions. 


I. INTRODUCTION 

Natural disaster^ such as floods or storms, represent 
extreme events with severe consequences including 
numerous killed and affected people as well as huge eco¬ 
nomic damage. Independent of the problem on how to 
project the occurrence of future extreme events, one is 
interested in which damage can typically be expected 
from an extreme event of certain magnitude. One ap¬ 
proach to tackle this question is to separate the statistics 
of extreme events from the damage caused by them. The 
former can be obtained from measurements, such as wa¬ 
ter level records, but the latter needs to be estimated by 
some empirical studies. 

Accordingly, so called damage functions provide a 
monetary value as a function of the magnitude of an 
event, such as the maximum flood level. We distinguish 
damage functions on a microscopic scale from those on 
a macroscopic scale The former describes the typ¬ 

ical costs of damages to single assets, such as residen¬ 
tial buildings, and the latter describes damage costs at 
larger areas, such as an entire city. This macroscopic 
(aggregated) damage function represents a composition 
of information on asset values, their location, and their 
vulnerability. We assume that in this case the damages 
are well characterized by the function, i.e. noise is re¬ 
duced due to spatial aggregation. 

In general, the functional form of damage functions is 
unknown. By definition, they are monotonic increasing 
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with the magnitude of the natural hazard and eventually 
exhibit saturation. For floods, on the microscopic scale 
a considerable set of recorded direct monetary damage 
values of inundated buildings has been related to the cor¬ 
responding water depth leading to an exponential depen¬ 
dence for private housing [ 1 , Fig. 5]. Furthermore, a set 
of different microscopic damage functions, models, and 
the related damages have been compared and evaluated 
Q , among them linear, quadratic, and square-root dam¬ 
age functions (see also [1^). In [Tl| a storm surge damage 
function has been estimated for the city of Copenhagen. 
In what follows we refer to macroscopic damage func¬ 
tions, i.e. the direct damage to an urban agglomeration 
as a function of the maximum flood height. 

Since from single sites there are usually too few dam¬ 
age records to obtain a damage function, we tackle the 
problem on a larger scale and study which functional 
form a damage function must follow so that the distribu¬ 
tion of extreme events transforms to the distribution of 
observed damages. Relating these two distributions we 
obtain macroscopic damage functions. Therefore, we an¬ 
alyze flood data assembled by CRED [I^ and find broad 
damage distributions. In order to characterize them, we 
elaborate two functional forms. We show that for the 
Gumbel case, Zipf’s law, i.e. a power-law distribution 
density with an exponent a 2, implies an exponen¬ 
tial damage function. Stretched exponential distribution 
densities imply power-law damage functions. 

As it is known analytically, maximum values of sam¬ 
ples of fixed size (such as ” block maxima” of time series) 
follow distributions which converge for sufficiently large 
samples towards Generalized Extreme Value (GEV) dis¬ 
tributions, if the values are independent and identically 
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distributed. Thus, we assume that changes in the statis¬ 
tics are small compared to the implied damages. GEV 
distributions are commonly fitted to maxima in order to 
estimate annualities and future occurrences [l^ . 

Further, if the data obeys GEV characteristics, it is 
possible to use the Generalized Pareto (GP) distributions 
as an approximation of the distribution function of the 
level s above a sufficiently high threshold st [a. There¬ 
fore, the upper tail of the distribution function of the 
(damage causing) events is described by one of the three 
GP distributions 


pOP^|l-(l+fi^) ' torl^O 
[l — e T for^ = 0. 

They are defined on s € [st, oo) and have a scale param¬ 
eter, 7 G K.+ , as well as a shape parameter, ^ € M (for 
^ < 0: s G {st,st — j))- According to the GEV-shape, 
one distinguishes three cases: (i) the Gumbel distribu¬ 
tion (^ = 0), (ii) the heavy-tailed Frechet distribution 
(^ > 0), and (hi) the bounded-tailed reversed Weibull 
distribution (^ < 0). For a more detailed presentation 
of extreme value assessment and applications we refer to 
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II. DAMAGE RECORDS 

We consider the EM-DAT database [l^ collected by 
the Gentre for Research on Epidemiology of Disasters 
(GRED) in version vI2.07 as created on Oct-28-2009. 
The information listed for each event entry consists of: 
start, end, country, location, type of disaster, sub-type, 
name, number of people killed, number of people affected, 
an estimated damage, and an ID. For the years 1950-2008 
we extract the information on floods, which include gen¬ 
eral floods, flash floods, as well as storm surges respec¬ 
tively coastal floods, and obtain 3469 entries worldwide, 
while for 1225 entries an estimated damage in units of 
Million US-Dollars is available. 

In Figure [iKa) we show the estimated probability den¬ 
sities, p{D), for all flood damage values of the database. 
However, since one may argue that the result could be 
biased by regional differences, in Fig. [IJb) we show p{D) 
for only those floods that occurred in the USA. Thus, we 
can weaken influences due to different economic power 
of different countries. In order to reduce possible trends 
in the data [l^, we also exclude floods before 1980 and 
plot the distribution density for Europe in Fig. [IJc). In 
any case, we observe broad distributions with damages 
reaching the order of 10 Billion US-Dollars. We would 
also like to note that we obtain similar distributions for 
the number of killed or affected people as well as for other 
natural disasters. 

In the simplest approach, the tail of the probability 
densities can be described with a functional form involv¬ 


FIG. 1. Probability densities of damages due to floods world¬ 
wide in the years 1950-2008 as estimated from the records pro¬ 
vided by CRED [^. (a-|-d) all entries, all years; (b-|-e) USA 
only, all years, 70 events; and (c-l-f) Europe only, 1980-2008, 
193 events. The straight lines in (a-c) are guides to the eye 
corresponding to Eq. m with a = 2. The solid lines in (d-f) 
are stretched exponential fits according to Eq. m providing 
the exponents (d) a ~ 0.32, (e) a ~ 0.42, and (f) a ~ 0.48. 
The probability densities have been estimated in logarithmic 
bins. 

ing one parameter, namely a power-law according to 

PiD) - , (2) 

where we find a « 2 [Fig. [U a-c)]. Such a size distribution 
is also known as Zipf’s law (a special case of the Pareto 
distribution) and is found in many different fields, such 
as word usage, city sizes, firm sizes, wealth, intensity 
of solar flares, etc. For an overview we refer to |19l - l^ 
and references therein. Minor deviations from Eq. @ 
for floods with small damage could be due to the fact 
that small damages are more likely to be missing in the 
database. 


III. RELATING EXTREME EVENTS AND 
DAMAGES 

The damage costs of a large flood magnitude event 
depend on a variety of damage influencing factors, such 
as orography, flow velocity, contamination, preparedness. 
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etc. [13 . Nevertheless, in the common practice, the dom¬ 
inant correlations with the maximum flood level are ex¬ 
plored [3]. Thus, if we assume a unique relation between 
the magnitude of an extreme event, s, and the mean dam¬ 
age, T*(s)j that is caused by it, we can take advantage of 
Eq. ([!]) and (l2|) and relate them. Therefore, we write the 
probability as an integral over the density and substitute 
the magnitude with the damage (s —>■ = Zl(s)): 

S2 ^(‘ 2 ) ^(‘ 2 ) 

Jp{s)ds= j J P(D)dD. (3) 

Si £>(si) £>(»i) 

Here, p(s) and p(£)) are the probability density of the ex¬ 
treme event and the damage, respectively. Furthermore, 

the density transformation = —— was used. Next, 

P{s) dD 

choosing si = st (the threshold for which GP distribu¬ 
tions are applicable) and S 2 = s, we obtain the equation 

D^s) 

P{s) = j P{D)dD, (4) 

which holds for any (reasonable) damage distribution 
density p(£)). 

Using Eq. @, i.e. the probability density P(d) = 
AD~°^ [a > 1) as indicated in Fig. [Ija-c), we write 

Pis)= J AD--dD = l-D<^-]Dl;^\ ( 5 ) 

^(st) 

where A = {a — 1)13“^“^ such that p(^D) is normal¬ 
ized in D G [D{st)^ P{s^oo))- Further we suppose 
that Il(s->.oo) 00 . Solving the equation for we 
get 

Pis) = D(st) (l “ Pis)) • (6) 

Finally, we insert the GP distribution for the Gumbel 
case, i.e. Eq. 0 with ^ = 0 , and obtain that the damage 
increases exponentially, 

Dis) = , ( 7 ) 

which holds for s G [stjOo). Accordingly, the exponen¬ 
tial damage function transforms the Gumbel distribution, 
approximated by an exponential one, into a Pareto dis¬ 
tribution [ 23 . 

However, the exponential damage function is based on 
damages following power-law distributions [Fig. [TJa-c)]. 
Since a one parameter description might not be sufficient, 
we also elaborate a two parameter fit - which we choose 
since it is integrable - namely a stretched exponential 
according to 

p{D) ^ , ( 8 ) 


where a and b are the parameters (a, 6 > 0). Equation dS]) 
is also known as Weibull distribution, see e.g. and 

references therein. For the same data as before the fitted 
curves are shown in Fig. [TJd-f) providing values for the 
exponent a roughly between 1/3 and 1/2. The stretched 
exponential (Weibull) distribution, Eq. dH), is sometimes 
used as penultimate approximation for pre- asymptotic 
behavior of extreme value distributions [^. This 
suggests, that the usage of the stretched exponential dis¬ 
tribution could be justified by an insufficiently high dam¬ 
age threshold. 

Now, the integral relating the extreme events s and 
damages D, Eq. 0 , is over p{D) from Eq. dll), instead 
of Eq. d2I): 


D. 


(o) 


Pis) = 


B-D^-h--!rdD =1-Be- 


(») 


(9) 


D 


(^t) 


(sr) ^ . 

where B = e i> such that P(_d) is normalized in H G 
[i:>(s,j,),i:)(s_,.oo))- Solving for £)(«) we find 


D 


(^) = 


D^^) - b ln(l 



( 10 ) 


Finally, we insert the GP distribution for the Gumbel 
case, i.e. Eq. m with ^ = 0 , and obtain 

Pis) = (PisT) + ■ ( 11 ) 

Accordingly, the power-law damage function transforms 
the Gumbel distribution, approximated by an exponen¬ 
tial one, into a stretched exponential distribution [^ . 
For a between 1/3 and 1/2 we obtain the asymptotic 
power-law relation Pis) ~ or Pis) ~ respectively. 


IV. SUMMARY AND DISCUSSION 


In summary, we characterize distributions of recorded 
flood damages, argue that they are caused by extreme 
events, and employ density transformation to deduce 
damage functions relating both. For Gumbel distributed 
extreme events (^ = 0 ) we find asymptotically 


Pis) 


3 

for p(£)) ~ D~°‘ with a > 1 
“ for p(£)) ~ witho> 0 . 


( 12 ) 

The involved GP parameters are in an ultimate sense and 
in practice penultimate approximations might be neces¬ 
sary [23[13|- In particular, it needs to be considered that 
if the underlying geophysical variable has only an approx¬ 
imate exponential distribution [i.e. ^ = 0 in Eq. O], then 
the form of the obtained damage functions would not 
be unique. An analogous argument applies to the cases 
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FIG. 2. Illustration of the obtained damage functions, (a) In the case of power-law distributed damages for ^ = 0 we 
obtain an exponential damage function [Eq. 0 , solid line]. For ^ 7 ^ 0 we insert Eq. ([T]) in Eq. ([ 6 )l and find faster (dotted) or 
slower (dashed) than exponential damage functions, Weibull (^ < 0) or Frechet (^ > 0) case, respectively, (b) In the case of 
stretched exponential distributed damages for ^ = 0 we obtain a power-law damage function with exponent 1/a [Eq. (Hill , solid 
line]. For ^ 7 ^ 0 we insert Eq. ([TJ in Eq. (IIOII and find that the damage function increases faster (dotted) or slower (dashed), 
Weibull (^ < 0) or Frechet (^ > 0) case, respectively, than the corresponding power-law. To generate the curves, we choose: 
Frechet, ^ = 0.1 and Weibull, ^ = —0.1 as well as a = 2 in (a) and a = 0.4 in (b). 


^ 7 ^ 0 and further research is necessary to find unique 
relations. 

The functional forms of Eq. (HU are illustrated in 
Fig. 0 which also includes the cases ■C 7 ^ 0- For power- 
law distributed damages and Weibull distributed extreme 
events (^ < 0 ) in average the damage increases faster 
than exponentially with the magnitude. For Frechet 
distributed extreme events (^ > 0 ) the opposite is the 
case. Intuitively, since in the Weibull case the extreme 
events have a less heavy tail, a steeper damage function is 
needed to result in the same damage distribution as the 
Gumbel case [Fig. [5][a)]. In the Frechet case, which has 
a fatter tail than the Gumbel distribution, a less steep 
damage function is sufficient to result in the same dam¬ 
age distribution as the Gumbel case. Similar arguments 
hold for stretched exponential damages [Fig. [21[b)]. 

We would like to remark that, allowing also negative 
values of a in Eq. ([5]), it represents the Frechet instead of 
the Weibull distribution. Moreover, with a scale param¬ 
eter a = — 1, in the limit s —>■ 00 the Frechet distribution 
is equivalent to Zipf’s law, Eq. m- 

In particular the role of aggregation on both, the ex¬ 
tremes and the damages functions requires further inves¬ 
tigations. While floods involve an integration of precipi¬ 
tation over a water basin, damages involve an integration 
of flood impacts over affected assets. Thus, it would be 
interesting to understand, to which extend the tail be¬ 
comes heavier at each stage, from precipitation extremes 
to flood magnitude and then to damages. 


Since the obtained damage functions represent an av¬ 
erage for large temporal and spatial scale, neglecting 
any further local differences, they have limited predictive 
power. However, the results could provide qualitative 
insight because exponential increasing damage is much 
more catastrophic than polynomial. Thus, it would be of 
interest to systematically analyze how the results depend 
on the considered spatial scale. 

Our approach could help to address the questions 
raised by decision makers or insurances, which costs cer¬ 
tain regions, countries, or even the globalized economy 
are facing from coastal floods [291 - K^ . The assessment 
of damage costs is fundamental in the econtext of adap¬ 
tation to climate change [s^. Finally, we would like to 
note that the presented approach can be applied to any 
type of natural disasters as long as the concept of damage 
functions can be used (see e.g. [H, H^)- 
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